%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%
% File:               predicted_probs.m
%
% Authors:            Sergio Ascencio and Miguel Rueda
%
% Language:           MATLAB R2013b (8.2.0.701) 64 Bit

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%


function [p_PAN,p_PRI]=predicted_probs(beta,X_PRI,X_PAN)

betas=reshape(beta,[size(X_PAN,2),4]);

%Base outcome is L (minimum representation) order of coefficients L, M, H
betas_PRI=[zeros(size(X_PRI,2),1) betas(:,1:2)];
betas_PAN=[zeros(size(X_PAN,2),1) betas(:,3:4)];

p_PAN=exp(X_PAN*betas_PAN)./kron(ones(1,3),sum(exp(X_PAN*betas_PAN),2));
p_PRI=exp(X_PRI*betas_PRI)./kron(ones(1,3),sum(exp(X_PRI*betas_PRI),2));

end

